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I. INTRODUCTION 

The development and analysis of continuum field theo- 
ries to model complex non-equilibrium structures or spa- 
tial patterns has made a tremendous impact in many 
areas of research in condensed matter and materials 
physics. A central idea in the development of such mod- 
els is the recognition that the patterns are controlled by 
the type and interaction of defects that define the pat- 
terns. For example in spinodal decomposition surfaces 
between different atomic species interact through diffu- 
sion in the bulk and control the phase segregation pro- 
cess. In block copolymer systems [l[ disclinations in- 
teract through elastic fields and control the ordering of 
lamellar phases. While it is desirable that models of these 
processes be derived from some fundamental atomic the- 
ory, they are frequently phenomenologically proposed. 
Classic phenomenological models include the Ginzburg- 
Landau theory of superconductivity and the Cahn- 
Hilliard-Cook theory of phase segregation 

Several years ago a phase field model of crystalliza- 
tion was phenomenologically proposed by exploiting the 
properties of free energy functionals that are minimized 
by periodic fields. In crystallization this field is inter- 
preted as the atomic number density (p), which is uni- 
form in a liquid phase and is typically periodic in the 
solid phase. By incorporating elasticity, dislocations and 
multiple crystal orientations, such functionals naturally 
incorporate the type and interaction of the defects that 
control many crystallization phenomena. This so-called 
phase field crystal (PFC) model 0, Q has been used 
to study glass formation Q, climb and glide dynam- 
ics [|[ , pre-melting at grain boundaries [H, [l(| , epitaxial 
growt h [111 Il2l . commensurate/incommensurate transi- 
tions [l3l. Il4|. sliding friction phenomena (l5j and the 
yield strength of polycrystals || [1, [HI]. More recently 
a simple binary phase field crystal modeled was devel- 



oped [T3| that couples the features of the PFC model of a 
pure materials with a concentration field so that eutectic 
growth, spinodal decomposition and dendritic growth can 
also be studied. This model can be linked with classical 
density functional theory (CDFT) (and the parameters 
entering CDFT), although the approximations needed to 
connect PFC with CDFT are quite drastic. As shown 
in recent studies on Fe [l8[ and colloidal systems [l^ |. 
CDFT predicts that p is very sharply peaked in space 
(at atomic lattice positions) while the PFC solutions are 
almost sinusoidal in space. Nevertheless these same stud- 
ies indicate that the parameters entering PFC models can 
be adjusted to match experimental quantities. 

While the periodic structure of PFC models is essential 
for describing elasticity and plasticity, it is very inconve- 
nient for numerical calculations. For example PFC simu- 
lations typically require 8 d (where d is dimension) spatial 
grid points per atomic lattice site. Obviously this limits 
the method to relatively small systems, although several 
new computational algorithms have been developed that 
can sig nificantly extend the applicability of both pure 
US El[ HH and binary PFC models To alleviate 
this limitation an amplitude expansion of the PFC model 
was developed by Goldenfeld et al. [25|, |2g, [27| . To un- 
derstand the idea behind such expansions it is useful to 
consider a one-dimensional equilibrium state of the form 
n = Asm(qx), where the amplitude, A, is zero in the liq- 
uid and finite in the solid state. While the field n varies 
rapidly in space, on a length scale set by A = 2ir/q, the 
amplitude A is uniform in crystalline regions and only 
varies near dislocations and liquid solid surfaces. De- 
formations of the crystal lattice can be represented by 
spatial variations in the phase of the amplitude. Using 
this amplitude representation, Athreya et al. [28| were 
able to apply adaptive mesh refinement to simulate grain 
growth on micron scales while simultaneously resolving 
atomic scale structures at interfaces. This remarkable 
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achievement suggests that the development of amplitude 
expansions is very promising for computational materials 
research. More recently this expansion has been extended 
to include spatial variations in the average number den- 
sity in two and three dimensions [29] ]. 

In addition to greatly increasing computational effi- 
ciency, amplitude representations of PFC models can also 
be exploited to establish a link between PFC type mod- 
els and traditional phase field models. This link provides 
insight into the specific terms that enter the bulk free en- 
ergy and gradient energy coefficients of traditional phase 
field models [3(], HH, [H, HH . Since the relationship be- 
tween the parameters that enter phase field models and 
sharp interface models are well established [34|, |H, [3(| , 
the relationship between parameters in PFC models to 
sharp interface models can then be established. 

In this paper, amplitude expansions are developed for 
triangular (2d), BCC and FCC crystal symmetries. The 
method of multiple scales expansion methods employed 
by Yeon et al. 29] is used. In the small deformation 
limit the expansions are shown to reproduce standard 
phase field models of solidification and eutectic growth 
which incorporate elasticity and solute segregation ef- 
fects. Sample simulations of grain boundary segregation, 
eutectic solidification and quantum dot growth on nano- 
membranes are also presented to illustrate the flexibility 
of the amplitude equations. 



II. BINARY PHASE FIELD CRYSTAL MODEL 

As discussed above, the binary-alloy PFC model devel- 
oped recently fl7j ] can incorporate the important features 
of solidification, phase segregation and solute expansion 
in alloy systems, in addition to the elasticity, plastic- 
ity, and multiple crystal orientations that characterize 
the crystalline state. For an alloy consisting of A and 
B atoms the model can be written in the case of equal 
atomic mobilities of the constituents as 
On 
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(M A + M B )/p 2 , M A (M B ) is mobility of 



where M = 

the A (B) species and r is time. The field n is the 
dimensionless number density difference given by n = 
(pA + Pb ~ Pe)/pe, where pa, Pb and pi are the atomic 
number densities of A atoms, B atoms and a reference 
liquid respectively. The field ip « (pa~ Pb)/ 'pe, plays the 
role of a concentration field. The free energy is given by 
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where A = 2i? 2 V 2 + i? 4 V 4 and t,v,~f,ui,u and K' are 
constants. The parameters B e , B x and R depend on i]j 



and in the simplest non-trivial case can be set to B l = 
B e + B^ 2 , B x = B x and R = R (l + a^). Details of the 
parameters entering this model are described in Ref. [13] • 
Here a is the solute expansion coefficient and to further 
simply calculations it will be assumed to be small. In the 
small a limit the free energy can be rewritten as 
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AB = B*-B%, F = AT/(k B Tp e R d a ), x = f/R a and 
K = K' /R 2 ,. The equations of motion are then, 
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where a time scale MksT piR d Q has been adopted for the 
rescaling. 

In the following section the equations of motion will be 
developed for the slowly varying amplitudes that describe 
various crystalline systems. More formally, for any given 
periodic structure the density field can be expanded as 
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which separates the "slow" -scale complex amplitudes rjj 
and average density field no from the underlying small- 
scale crystalline structure that is characterized by G — 
Iqi + m<j2 + nqz, where ((71,92,93) are the principle re- 
ciprocal lattice vectors and (l,m,n) are integers. The 
corresponding solutions for n are then relatively smooth 
and can be approximated using only a few amplitudes. 
In what follows, model equations will be developed for 
only the lowest order amplitudes that are needed to re- 
construct a given crystal symmetry and defect structures 
relevant for controlling elastic and plastic effects in solid- 
ification and impurity segregation. 

Recently many efforts have been devoted to developing 
amplitude expansion for various physical systems. The 
central assumption of these approaches is that the am- 
plitudes vary on scales much larger that the short (or 
"fast") atomic spacing scale. Formally a small parame- 
ter can then be introduced that represents the ratio of 
the two scales and an expansion in this variable can be 
performed. To apply this analysis to Eqs. ([5]) and ©, 
both the amplitudes (77^- ) and concentration field ip are 
assumed to be slow variables. For technical details of 
this multiple-scale analysis the reader can refer to Yeon 
et al. [29l | and the references therein. For simplicity in 
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this paper the average atomic density (i.e., no) will be 
assumed as constant and zero, since miscibility gaps be- 
tween liquid and solid phase can be accounted for (to 
some extent) by a miscibility in ip. Moreover, noise dy- 
namics will be neglected here. More complete analysis 
involving dynamic variation of no, noise effects, as well 
as the general case of different atomic mobilities will be 
presented elsewhere [37| . 



III. AMPLITUDE EXPANSION FOR 
TRIANGULAR SYMMETRY 

In two dimensions the free energy given in Eq. ([3|) 
is minimized by a triangular lattice. The corresponding 
principle reciprocal lattice vectors are given by 



Qi = qeq 



(-x - 1/V3yj , q 2 = q eq (2/V3yj , (8) 



limit the equations of motion become 
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These dynamics can alternatively be written in a form of 
Model C ty pe (in the Halperin-Hohenberg classification 
scheme |3{|), i.e., 



where q eq — \[Zj2 is the equilibrium wave number. To 
construct the minimal model of a triangular lattice only 
the lowest order reciprocal lattice vectors are needed, 
which correspond to gi, q\ and q\ = —q\ — q\ = q e q{x — 
l/v3$)- Follo wing the standard methods of multiple- 
scale expansion [29|, [38| the following equations of motion 
can be derived 
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where 



with c.c. referring to the complex conjugate. 
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A) = ABgCjQj 

Cj = V 2 + 2iqj ■ V - \q 3 \ 2 
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3 = 1 

To further simply calculations a long wavelength approxi- 
mation will be made such that C j « — | qj \ 2 = — 1 . Unfor- 
tunately a similar long wavelength approximation can't 
be made for Gj (e.g., replacing Gj by 2icfj ■ V) as then the 
equations would not be rotationally invariant. In this 



A. Small deformation limit 

To gain insight into the above results and make con- 
nection with traditional phase field models it is useful 
to rewrite rjj in the form rjj = <fie iqj ' u . In this case 
the magnitude of (p distinguishes between liquid {(p = 0) 
and solid phases {<p ^ 0), while u is the displacement 
vector introduced in continuum elasticity theory [42| to 
describe displacement of atoms from a perfectly ordered 
crystal lattice. Substituting this expression into Eq. (fT3")) 
and taking the long wavelength limit (done by retaining 
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derivatives only to second order) gives, 
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where Uij = (dui/dxj + duj/dxi)/2 is the linear strain 
tensor. The first set of terms (defined by the {...} brack- 
ets) is remarkably similar to standard phase field models 
developed for eutectic and dendritic solidification [40L[4ll j . 
The polynomial in <f> gives a first order transition from a 
liquid {4> = 0) to solid phase (0^0). The polynomial in 
-0 is the typical l- 4 ' free energy used in the Cahn-Hilliard- 
Cook type models. The coupling term (f> 2 ijj 2 (note: B\ is 
negative) can lead to phase segregation at low tempera- 
tures when <fi becomes large. The second and third set of 
terms account for surface and linear elastic energy respec- 
tively. From the form of the third term it is straightfor- 
ward [42j to derive the elastic constants in dimensionless 
units, i.e., C\\ = C22 = 9i?Q0 2 and C\i — C44 = Cn/3. 

Finally the last set of terms couples the concentra- 
tion field to the liquid/solid order parameter <p when 
the atomic species have a different size (i.e., a =/= 0). 
The term, ip(f>'V 2 (f>, implies preferential phase segregation 
to liquid/solid surface, grain boundaries and dislocations 
(i.e., regions in which <p varies in space). Dynamically this 
term is responsible for solute migration at grain bound- 
aries and solute drag. The last term ip(f> 2 (U xx + U yy ) 
implies a coupling between strain and concentration as 
should be expected when the atomic species have differ- 
ent sizes. It may appear unusual that the free energy 
functional depends on the sign of the displacement gra- 
dients (via U xx = du x /dx); however this sign determines 
whether there is a local compression or expansion of the 
lattice which would favor a specific atomic species based 
on the sign of solute expansion coefficient a. 



o 
o 



CQ 

-a 




FIG. 1: Sample phase diagram for two dimensional tri- 
angular system with parameter set (Bq , B|, v, t, a, it) = 
(1, -1.8, 1, 3/5, 0.30, 4). Also lo = 0.088 and 0.008 in (a) and 
(b) respectively. In each panel the filled regions correspond 
to regions of phase coexistence. 



Minimizing with respect to 8 gives S eq = —aip. As ex- 
pected the contraction/expansion of the lattice is con- 
trolled by atjj. Substituting S eq for S leads to 
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which is then minimized with respect to <fi, yielding 
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Substituting Eq. (Tl~7|) into Eq. (fTf]]) yields a free energy 
per unit area that is only a function of ip. This result can 
be used to construct the phase diagram in an analogous 
manner as was done in Ref . (l7| . Two example phase dia- 
grams and the corresponding model parameters required 
to obtain each are show in Fig. [TJ 

The derivation of amplitude equations and the related 
simplifications presented above can be readily extended 
to three dimensional systems. In Appendix [Al and [Bl the 
relevant equations and sample phase diagrams for both 
BCC and FCC symmetries are presented. 



IV. APPLICATIONS 



B. Equilibrium phase diagram 

The equilibrium phase diagram can be evaluated by 
considering <f> and ip constant and a bulk compression to 
account for solute expansion, such that u = d(xx + yy). 
In this limit the free energy per unit area (A) becomes 

I = 3AB o 2 -4^ 3 + ^ 4 + ^ 4 

,/,2 

+ (w + 6B e 2 </) 2 ) ^- + 12B*ct) 2 (5 2 + 2aip6) . (15) 



To further examine the above amplitude equations, nu- 
merical simulations of the model described by Eqs. (jX3J) 
and (|13|) were undertaken and the results are shown in 
Figs. [2HSl In Fig. [2] simulations were conducted to study 
the coupling between composition and topological defects 
(dislocations and grain boundaries) in binary alloys with 
components of different atomic sizes, i.e., nonzero solute 
expansion coefficient a. For this study a symmetric tilt 
grain boundary between two grains with a misorientation 
angle of 9 = 3.76° was created by dynamically evolving 
an initial configuration of two perfect crystals separated 
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FIG. 2: Solute segregation in symmetric grain boundaries 
with misorientation 8 = 3.76° are shown for tp = (a) and 
■0 = 0.2 (b). The left and right panels correspond to ■ |»7j| 2 
and tjj respectively. In the corner insets of left panels the 
density field n is reconstructed from the amplitudes rjj for 
the boxed region. In the right panels the dark (light) color 
corresponds to the larger (smaller) of the atomic species. The 
parameters for (a) and (b) correspond to Fig. [T^, and Fig. [lb 
at ABq = 0.01 respectively. 



by a layer of supercooled liquid. As time evolves the liq- 
uid solidifies and a grain boundary spontaneously forms. 
The dislocation cores that comprise the grain boundary 
interact with the different atomic species or alloy com- 
position. As shown in Fig. [2] the larger (smaller) so- 
lute atoms preferentially segregates around the disloca- 
tion cores in regions of tensile (compressive) strain. 

Simulation results for eutectic solidification and phase 
separation are presented in Fig. [3J where three small 
grains of different orientation are heterogeneously nucle- 
ated for the parameters used in Fig. [TJd and for a "tem- 
perature" AB below the eutectic point. As the grains 
grow the lamellar concentration bands form within the 
grains as a result of phase separation. The relatively large 
lattice mismatch (roughly 8.4% in equilibrium, due to a 
finite solute expansion coefficient) between the lamella 
results in the spontaneous nucleation of dislocation at 
the lamellar interfaces. When the grains impinge on one 
another (i.e., coarsening occurs) additional dislocations 
form, leading to complex patterns as shown in Figs. \Sjp 
and[3j;. All these simulations indicate that the amplitude 
model established here can simultaneously describe the 
complex evolution of liquid/solid interfaces, grain bound- 
aries, dislocations and interfaces between regions of dif- 
ferent compositions. 

Numerical simulations were also conducted to study 



FIG. 3: Eutectic solidification for parameters in Fig. [Tb at 
ABo = 0.022 and ip = 0.0. Panels (a), (b) and (c) correspond 
to dimensionless times 30000, 60000 and 105000 respectively. 
From left to right the columns correspond to n reconstructed 
from Tjj (for the boxed region), . |% | 2 , 4> an d the local free 
energy density. Dislocations are most easily identified as small 
black dots in the local free energy density. 



islands or quantum dots formation on thin freely stand- 
ing films or solid nanomembranes. Recent experiments of 
Ge on Si membranes [H, |44[ have suggested that growth 
on such nanomembranes strongly influences the maxi- 
mum size that the strained islands can form coherently 
(i.e., without dislocations) and can lead to correlation or 
self-ordering of multiple islands. To examine this phe- 
nomenon, simulations were set up such that islands of 
one material were grown on a thin free standing mem- 
brane of another material, by exploiting a eutectic phase 
diagram (such as that shown in Fig. QJ>). To initiate 
growth of the islands a small crystal region at i/j = —0.12 
was constructed coherently on top of a thin membrane at 
-0 = 0.2 in a supercooled liquid at ip = —0.03, with AB 
set to be 0.028. Sample simulation results are shown in 
Figs. 2] and [5] As can be seen in these figures the dot 
grows coherently with the membrane until reaching some 
critical size at which dislocations are nucleated at the liq- 
uid/dot/membrane junction. Perhaps more interestingly 
by comparing Figs. 0] and [5] it is apparent that for the 
thinner membrane the dots can grow to a larger size be- 
fore nucleating dislocations. This result occurs as the 
strained quantum dots partially relax by straining the 
substrate membrane and thinner membranes are easier 
to deform than thicker ones. Such mechanism has been 
proved to play an important role in engineering the self 
assembly of thin film nanostructures such as quantum 

dots mf. 

Another interesting consequence of the membrane de- 
formation is that it locally creates favorable and unfa- 
vorable positions for the nucleation of other islands/dots, 
with an alternative sequence on the two sides of the mem- 
brane, as can be seen in Fig. [5] which was conducted at a 
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FIG. 4: Quantum dot formation on a two atomic layer thick 
nanomembrane at times t = 20000, 60000 and 100000 (for 
(a) to (c) respectively). The columns from left to right cor- 
respond to ^2 \Vj\ 2 j V' an d the local free energy density. The 
parameters used for this simulation are from Fig. [ij, except 
a = 0.26 and B l = 1.028. 
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FIG. 5: Quantum dot formation on a five atomic layer thick 
nanomembrane at times t = 20000, 40000 and 60000 (for (a) 
to (c) respectively). The parameters used are identical to 
those in Fig. |4] 




FIG. 6: Correlated quantum dot formation on a nanomem- 
brane at times t = 20000, 40000, 60000, 80000, 100000 and 
120000 (for (a) to (f) respectively). The parameters used 
are identical to those in Fig. U except for a higher liquid 
supersaturation. 



higher liquid supersaturation (i.e., ip = —0.04) or growth 
rate. Once a dot is formed, e.g., on the top of the mem- 
brane, it is preferable for the next ones to nucleate under 
the edges rather than directly underneath the top quan- 
tum dot. After these bottom-side dots are formed they 
in turn create preferential regions for a new set of dots to 
nucleate on the top (above their edges) and the process 
repeats. In such a fashion the process leads to the cor- 
related dot growth as observed in experiments of Ge or 
SiGe islands on Si membranes [43, 44]. This phenomenon 
can in principle be exploited to create periodic strained 
nanostructures that in turn produce periodically modu- 
lated band gaps. 

Finally it should be noted that in the dynamics de- 
scribed by Eq. (fH?|) the diffusion constant for the con- 
centration field (-0) is similar in magnitude in the liquid 
and solid phases. However it is typically the case that dif- 
fusion or atomic mobility in liquids is much larger than 
in solids. This difference can be incorporated by using 
the following dynamic equation for ip, 

where the mobility now depends on A 2 = J^ . \r)j\ . The 
following functional form, 

M = M x + (M L - M x ) [1 - tanh (aA)] (19) 

would for example change M from the mobility of the 
crystal (Mx) at large aA to the mobility in the liquid 
phase (Ml) at small or zero aA. The parameter a would 
control how quickly the M changes from liquid to solid 
phases. In addition, while the binary PFC model pre- 
sented in Eq. fl} assumes both atomic species have the 
same mobility, the concentration dependence of M could 
also be added in the above ad hoc method suggested for 
the phase dependence of atomic mobility. 

V. DISCUSSIONS AND CONCLUSIONS 

In this paper amplitude equations have been derived 
from the binary phase field crystal model for a two di- 
mensional triangular lattice and for three dimensional 
BCC and FCC structures. Furthermore, the connec- 
tion to standard phase field models has been established, 
and for small deformations the results have been shown 
to recover linear continuum elasticity theory and recon- 
struct the equilibrium phase diagrams for binary alloy 
systems. Sample simulations of the amplitude equations 
have shown that this relatively simple model can effec- 
tively model many complex phenomena and the emer- 
gent microstructures that arise, and reveal the underly- 
ing mechanisms. While these amplitude equations were 
derived from an atomistic model, they can in themselves 
be regarded as phase field models with complex order pa- 
rameters. One advantage of this amplitude description is 
that the liquid and solid phases are easily distinguished 
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by a relative uniform quantity A 2 = ^ \rjj\ 2 , and the 
coupling and interaction between structure (i.e., ampli- 
tudes) and concentration can be well identified. In this 
regard it may be possible to extend the equations to natu- 
rally incorporate other uniform fields, such as magnetiza- 
tion, polarization, temperature, etc. and simultaneously 
include elastic and plastic deformation in polycrystalline 
samples. 
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APPENDIX A: AMPLITUDE EQUATIONS FOR 
BCC SYMMETRY 

For a BCC lattice, the principle reciprocal lattice vec- 
tors are 



qi = q eq (x + y), 
eh = q eq (x + z), 
93 = q eq (y + z), 



(Al) 



where q eq — l/y/2 is the equilibrium wave number. In a 
one mode approximation the above principle reciprocal 
lattice vectors need to be combined with the following 
vectors, q 4 = q\-q\ = q eq (y-z), $5 = <72-<?3 = Qeqix-y) 
and <f 6 = <f 3 — q\ = q eq (—x + z). Thus the BCC structure 
can be represented in the usual manner, i.e., 
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Repeating the calculations presented in section (jllll 
(with the same level of approximations) gives the fol 
lowing complex amplitude equations: 
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with equations of motion for 772 and 773 obtained by cyclic 
permutations on the groups (1,2,3) and (4,5,6) in Eq. 
[3l) . while equations for 775 and r]Q can be obtained by 



FIG. 7: Sample phase diagram for three dimensional BCC 
system with same parameter set and notation as in Fig. [1] 
except u> = 0.05 in (a). 



the similar cyclic permutations of Eq. (|A4[) . The corre- 
sponding concentration equation is given by 



~dt 



V 2 ((w + B e 2 A 2 - KV 2 ) ip + u^ 3 



(A5) 



Once again, Eqs. (|A3 
'model C type form, i.e., 



5| can be rewritten in a 



dt 



SF 



~dt 



,5F 



(A6) 



where 



F 



Ago ^2 



3i> 



A 4 



El^lW-yl^l 4 } 
3=1 L J 



+6tj (mrisVlvl + + m^VeVt + cc.) 

-2t ([rjlmVi + V2V3V5 + V3V1V6 + cc] 
+ [vtVsVe +c.c.]) 



3=1 



(A7) 



As discussed in Sec. MI Al it is interesting to replace rjj 



8 



by §e lqyr , yielding 



F = J dr[{6AB Q (j) 2 -16t<j) 3 + I35v<j) 4 



■0 2 



[i=i \ j¥» / 

6 1 

+ 2 E^ 



(A8) 



i=4 



Similarly the elastic constants can be derived as Cn = 



-13 



G 



2:; 



Caa — 



C22 = C33 = 8B'§^ 2 and C 12 = C x 
C55 = 6*66 = Cii/2. Furthermore, minimizing with re- 
spect to S gives, S eq = leads to the following result 
for the free energy per unit volume, 



- = 6AB O 2 - 16i0 3 + 135u0 4 



+ {lo + &B l 2 4> 2 ) 



w 2 



+~0 4 -24B o *(c^>) 2 ) 



(A9) 



which can be minimized with respect to (f> to obtain 
_ 2t + y/4t 2 - A5v(AB + ijj 2 (B e 2 - 4B x a 2 ) 

= A 10) 

The corresponding phase diagram for this model is shown 
in Fig. [3 



APPENDIX B: AMPLITUDE EQUATIONS FOR 
FCC SYMMETRY 



Recently Wu [45| has shown that the basic phase field 
crystal model can be extended to model a FCC lattice by 
including an extra length scale. Replacing the operator 
AB + B x (1 + i? 2 V 2 ) 2 -"►£> = AB + 16B X (1 + 2R 2 1 V 2 + 

i? 4 V 4 )(l + 2i? 2 V 2 + i? 4 + £0, where R 2 = \/3/2Ri for an 
FCC lattice and the extra factor of 16 is introduced just 
for convenience, gives a revised alloy PFC free energy, 



F = 



dr 



riDn t 



n 

3 



— n 

4 



(Bl) 



Here the parameter E controls the symmetry of the phase 
such that for large (small) E a BCC (FCC) structure is 
favored. The following calculations focus on the FCC 
phase in the limit of E — 0. Furthermore, to examine 
the influence of solute expansion the parameters R\ and 



i?2 can be set to be R\ = 1 + aip and R 2 = \/3/2(l+a')jj). 
In the small aijj limit the free energy functional becomes, 



F = dr 



lA°7 



aip- 



—n 
3 



— n 

4 



K 



(B2) 



showing the same form as Eq. © but with different 
operators 

A° ee AB + B^ 2 + 16B X {1 + V 2 ) 2 (1 + 'S/4V 2 ) 2 , 



A 



\6B X V 2 {\ + V 2 )(l + 3/4V 2 )(7 + 6V 2 ). (B3) 



The corresponding equations of motion are also governed 
by Eqs. ([5]) and ©, with A and A 1 given above. 

For the FCC symmetry the principle reciprocal lattice 
vectors are 



Ql = Qeg(-X + y + £)/v / 3, 

92 = q eq {x - y + z)/Vs, 
<?3 = qeq{£ + y- z)/Vs. 



(B4) 



To construct an FCC crystal the following reciprocal lat- 
tice vectors are also required, 

9*4 = -<?i-<?2-<?3 = qeq(-z-y-£)/V3, 

95 = Ql + ?2 = 2q eq z/V3, 

96 = 92 + 93 = 2q eq x/ \/3, 

97 = 93 + 9i = 29eg2// V3. (B5) 

Unlike the triangular and BCC symmetries, the FCC lat- 
tice requires at minimum two set of vectors of different 
lengths, i.e., (q x , <f 2 , 93, 94) with length q eq and (95, q 6 , 97) 
with length 2/y/3q eq . The density field n is then ex- 
panded in the usual fashion, i.e., 



3=7 
3=1 



(B6) 



Following the standard procedure, the amplitude equa- 
tions can be derived as 



drji 

at 



d [A 1? /i - 2t (rj*r] 5 + 77*777 + 7]* 4 r]*) 
+6v({A 2 -\r ll \ 2 /2}r ll +r 1 * V * 3 r ] * 4 



l 2 ~lm| 2 

+V2V6V7 + V3V5V6 + %?75?77) 



a 
2" 



(B7) 



dr] 2 

at 



£2 [A° 2 ?/2 - 2t (^3776 + r)Xrfr + r\\r) h ) 
+6v({A 2 -\ m \ 2 /2} m +r 1 * lV * V l 
+V1V6V7 + 773775777 + 774775776) 

+ 77 



(B8) 
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o 

2 2 

x 



m 

< 1 



-v v i 1 1 1 1 i 1 1 1 1 y v- 

- \.\. liquid / y '_ 
solid 


-1 'MM 1 1 1 1 1 111 yv - 
- \ \ liquid / / '_ 

' SOlld^ff^;)''':^, M)M : ':k Solid " 




1 li:i!ii:i!ii:i!iii:i!ii:iniUni:inV:ini::i!k 



d^_ 
~dt 



V 2 ((w + B\A 2 - KV 2 )ip + uijf 



a 



Y^bA^YV'+V-^M , (B14) 



where C 3 = V 2 + 2iq 3 ■ V - \q 3 \ 2 , A 2 = 2 V ■ \Vj\ 2 , and 



-0.05 0.05 -0.1 -0.05 0.05 0.1 

FIG. 8: Sample phase diagram for three dimensional FCC 
system with the same parameter set and notation as in Fig. 
Q] except lu = 0.02 in (a). 



dm 

dt 



dm 

dt 



dm 
dt 



dm 
dt 



dm 
dt 



£ 3 [A° 3 773 - 2t (ritvt + V*m + V* 2 Ve) 

+6v ({A 2 - \m\ 2 /2} m + vMvl 
+vivlve + mnim + mmm) 



(B9) 



U [A 4 r74 - 2t (vlVe + V* 2 V*7 + vH) 

+6v({a 2 -\m\ 2 /2}m + v*iv* 2 m 
+mvlv7 + V2V5V6 + V3V6V7) 

+ -(V'A 1 4 ?74 + A 1 4 774) 



£5 [A5775 - 2t {mm + V3V4) 

+6v ({A 2 - \m\ 2 /2} m + vemvs 
+mmm + mmm + mmm) 

(^A\m + ^\m)} , 



a 
2 



(Bll) 



£ 6 [A° 6 7 76 - 2t (mm + vtvl) 
+6v ({A 2 - \m\ 2 /2} m + mmvi 
+V7mvl + mmvl + vlmvl) 

(V'Aerye + A^e)] , 



a 
2 



(B12) 



c 7 [A° 7 m~2t (mv3 + V2vl) 

+6v({a 2 -\m\ 2 /2}m + V5v*2V3 
+V5V1VI + mmm + mmm) 



a 



(i>k\m + aW) 



(B13) 



A) = AB + B l 2 ^ 2 + 16B X {1 + C 3 ) 2 (l + 3£ 3 /4) 2 , 

16B X £ 3 (1 + Cj)(l + 3/4£j)(7 + 6C 3 ). (B15) 



4 



Note that \q 3 \ 2 = 1 for j = 1,2,3,4 and |q.,-| 2 = 4/3 for 
j = 5, 6, 7, and thus we have 



A° = 



AB + BkU? + /Mc/ ;i -(3c;, + l) 2 j = 1, 2, 3, 4 



AB + B^j 2 + B X (3G 3 - lf{Qjf j - 5,6, 7 
AB + B^ 2 + B x (g 3 ) 2 j = l,..., 7 (B16) 



A 1 ,- 



4 K (£y - 1)^(3^ + l)(6£y + 1) 3 = h2,3,4 

45^ - 4/3) (3^ - 1)^(6^ - 1) j = 5, 6, 7 

-4B*(0 3 -) i= 1,2,3,4 
-4B*(40i/3) j = 5,6,7 



(B17) 



where C?j = V 2 + 2ig} • V. Finally as described in Sec. 
IIII Al the operator Cj acting on the right hand side of 
each equation can be approximated by constants, i.e., 

(BIO) 

Cj ~ — \q 3 \ 2 ■ As before the dynamical equations can be 



written as 



dm _ _ w ,sf 
dt " m s v y 



dip 



,SF 



dt 



(B18) 



where 
F = 



dr 



-A* 



j^W\QM 2 ^\mt\ 

.7=1 L ' 



-2t (vl(v2V5 + vim + vtve) + vlinlvs + vlrf) 

+Vlvtv* 5 + c.c.) + 6v (VKV2V3V4 + V2V6V7 

+mm~nl + mmm) + iilmimrf + mm) 
+t]lmmm + c -c-) 



J =1 



+ ^ + BiA 2 )^ + ^-^\Y i - 



3=5 

K, 
~2 



(B19) 
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Setting r]j = 4>e lqj ' r as before gives 

693w 



F = Id? 



7AB Q( f) 2 - 2Ucj) 3 



~K~ B o { E I W l + Ua ^ + E U " U U 



-4E^ 

i=4 




(B20) 



from which the elastic constants for a FCC lattice can 
be obtained: C n = C 22 = C 33 = 160 B§<j> 2 / 9 and Ci 2 = 



C13 — C23 — C44 — C55 — — 32_Bq0 2 /9 — Cii/5. 

If minimizing the above free energy expression with re- 
spect to 5, we obtain 5 eq = —mp and hence the following 
free energy per unit volume 



- = 7AL'„o 2 -24/o : + ^ u ^ + ( w + 7fl|0 2 )^ 



F 
V 



+^ 4 -^BJ(a#) 2 , 



2 

(B21) 



which is minimized when 



18< + v /324i 2 - 4851w(AB + ip 2 (B l 2 - 16^a 2 /3)) 



693v 



(B22) 

A phase diagram for this FCC amplitude model is pre- 
sented in Fig. [5] 
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